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Electron - hole pairs are copuously created by an applied electric field near the Dirac point in 
graphene or similar 2D electronic systems. It was shown recently that for sufficiently large electric 
fields E and ballistic times the I-V characteristics become strongly nonlinear due to Schwinger's 
pair creation rate, proportional to E^^'^ . Since there is no energy gap the radiation from the pairs' 
annihil ation is enhanced. The spectrum of radiation is calculated and exhibits a maximum at 
u] — \J eEvg/h. The angular and polarization dependence of the emitted photons with respect to 
the graphene sheet is quite distinctive. For very large currents the recombination rate becomes so 
large that it leads to the second Ohmic regime due to radiation friction. 

PACS numbers: 72.80.Vp 73.20. Mf 12.20.-m 



INTRODUCTION 

Electronic mobility in graphene, especially one sus- 
pended on leads, is extremely large [l[ so that a graphene 
sheet is one of the purest electronic systems. The relax- 
ation time of charge carriers due to scattering off impuri- 
ties, phonons, ripplons, etc., in suspended graphene sam- 
ples of submicron length is so large that the transport is 
ballistic [H, Q . The ballistic flight time in these samples 
can be estimated as thai = L/vg , where Vg ~ lO^m/s is 
the graphene velocity characterizing the massless "ultra 
- relativistic" spectrum of graphene near Dirac points, 
£fe = Vg |k|, and L is the length of the sample that can 
exceed several /xm 0, IH . The extraordinary physics ap- 
pears right at the Dirac point at which the density of 
states vanishes. In particular, at this point graphene ex- 
hibits a quasi - Ohmic behaviour, J = crE, even in the 
purely ballistic regime. 

A physical picture of this "resistivity" without either 
charge carriers or dissipation is as follows Q . The electric 
field creates electron - hole excitations in the vicinity of 
the Dirac points similar to the Landau-Zener tunneling 
effect in narrow gap semiconductors or electron - positron 
pair creation in Quantum Electrodynamics first studied 
by Schwinger 0] (later referred to as LZS). Importantly, 
in graphene the energy gap is zero, thus the pair cre- 
ation is possible at zero temperature and arbitrary small 
E, even within linear response. Although the absolute 
value of the quasiparticle velocity Vg cannot be altered 
by the electric field due to the "ultra - relativistic" dis- 
persion relation, the orientation of the velocity can be 
influenced by the applied field. The electric current, ev, 
proportional to the projection of the velocity v onto the 
direction of the electric field is increased by the field. 
These two sources of current, namely creation of moving 
charges by the electric field (polarization) and their re- 



orientation (acceleration) are responsible for the creation 
of a stable current. 

Agreement over the qualitative explanation notwith- 
standing, determination of the value of the minimal DC 
conductivity at Dirac point in the limit of zero tempera- 
ture had undergone a period of experimental and theoret- 
ical uncertainty. After the value in graphene on substrate 
was measured to be about a — Ae'^/h [8|, it was shown in 
experiments on suspended samples that the zero tem- 
perature limit was not achieved and in fact that these 
early samples had too many charged " puddles" , so that 
they represented an average around the neutrality or the 
Dirac point. The value in early-on suspended samples 
was half of that and most recently settled at the "dy- 
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namical" (72 = f x hest samples at 2K temperature 
3 . Theoretically several different values appeared. The 
value 0"! = had been considered as the "standard" 
one for several years Q and appeared as a zero disorder 
limit in many calculations like the self consistent har- 
monic approximation, although different regularizations 
within the Kubo formalism resulted in different values 

The dynamical approach to transport was applied to 
the tight binding model of graphene [ll| to resolve this 
" regularization ambiguity" . It consists of considering the 
ballistic evolution of the current density in time after a 
sudden or gradual switching on of the electric field. The 
result within linear response is that the current settles 
very fast, on the microscopic time scale of = h/j ~ 
0.24 fs (7 being the hopping energy), on the value of 
J = a2E. The value is identical to the one obtained (at 



nonzero temperatures) for the AC conductivity [12[. The 



two contributions, polarization and attenuation are com- 
parable in strength and combine to produce a constant 
total current. However a deeper analysis of the "quasi 
- Ohmic" graphene system beyond the leading order in 
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perturbation theory in electric field revealed |l3| that on 
the time scale 



tnl 



eEvn 



(1) 



the linear response breaks down. For larger times the 
quasi - Ohmic behavior no longer holds. This is in con- 
trast to dissipative systems, in which the linear response 
limit can be taken directly at infinite time. This per- 
haps is the origin of the " regularization" ambiguities in 
graphene, since large time and small field limits are dif- 
ferent. The time scale on which nonlinear effects become 
dominant is not always very large; for example, in ex- 
periments dedicated to breakdown of Quantum Hall ef- 
fect [31 in which E = W^V/m, nonlinearity sets in at 
tnl = 0.3ps, that is of order ballistics time for L = 0.3/im. 
Graphene flakes under larger fields of order 2 • lO^V/m 
have been studied very recently (at room temperature) in 
specially designed high current density experiments Q. 
In this case the nonlinear time is only 20/s, much lower 
than the ballistic time tbai = ^ps for L = 2/im. Analytic 
and numerical solutions of the tight binding model[13l|. 
as well as of the Dirac model describing the physics near 
the Dirac point demonstrated 13, 1^ that at tni the elec- 
tron - hole pairs creation becomes dominant and is well 
described by an adaptation of the well - known (non- 
analytic in E) Schwinger electron - positron pair creation 
rate 



-Nr, 
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eE 
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The difference with the original derivation 0] in the con- 
text of particle physics is that the fermions are 2+1 di- 
mensional and "massless", thus magnifying the effect. 
The polarization current is J [t) ~ 2eVgN (t) and there- 
fore Schwinger's creation rate leads to a linear increase 
with timell3l: 



J (t) = a2 



■^/3 



3/2 



E 



1/2 



(3) 



The physics of pair creation is highly non-perturbative 
and non-linear in nature and therefore, instead of the 
linear response, Schwinger found an exact formula using 
functional methods. The rate can be intuitively under- 
stood using the much simpler instanton approach origi- 
nally proposed in the context of particle physics [3 (ex- 
tended later to low dimensions [l7| ) , but is known in fact 
in condensed matter physics as the Landau - Zener tun- 
neling probability 0, [l5, 18|. In particle physics it is ex- 
tremely difficult to observe Schwinger's creation rate and 
it would be interesting to establish experimentally this 
dynamical phase in low dimensional condensed matter 
physics featuring the massless Dirac quasiparticle spec- 
trum like graphene or novel materials sharing with it 



the massless Dirac spectrum like topological insulators 
or tuned semiconductor heterojunctions [l9|. Of course 
transport phenomena at rather large fields always have a 
background related to possible influence of leads, disor- 
der and thermal effects like local heating, etc. 

In this note we draw attention to a direct and unin- 
trusive signature of the dynamical phase of LZS pair cre- 
ation in a graphene sheet subject to an applied electric 
field. It is demonstrated that the flux of photons radi- 
ated from the surface of the sample is characterized by 
the creation rate since the photons are emitted via elec- 
tron - hole pair annihilation and therefore proportional 
to iJ"^/2, a hallmark of Schwinger's process. In addition, 
the frequency, direction and polarization characteristics 
of the radiation generated by the electric field calculated 
here all bear footprints of the pair creation dynamics. 



ELECTRON - HOLE RECOMBINATION RATE 
INTO PHOTONS. 

Amplitude for emission of a single photon 

The electrons and their electromagnetic interaction 
with photons are approximately described near a Dirac 
point by the Weyl Hamiltonian: 



H 



-ihV 



Vr_o 



(4) 

Here if) is the two component spinor second quantized 
field and (A,^^) is the vector potential (bold letters de- 
scribe vectors in the graphene plane, while z is the per- 
pendicular direction). Electrons (charge — e) and holes 
(charge e) in the graphene sheet are confined in this 
model to the 2: = plane by a potential Vconf (small 
shape changes can be neglected for our purposes). The 
only requirement from this potential is that it is strong 
enough to "freeze" the motion along the z direction. In 
the single graphene sheet one has two left handed chiral- 
ity Weyl fermions described by the above Hamiltonian 
in which cr denotes the in - plane Pauli matrices and 
two right handed Weyl fermions represented by cr'^. To 
include the topological insulators case [l9l |. we first con- 
cern ourselves with only one spinor. 
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Fig.l. Feynman diagrams representing the major 
electromagnetic processes in graphene. a. The one 
photon emission. 



p-k 





b. The two photon emission 



We consider the emission of a photon with wave vector 
(k,fcj;) and frequency cj = c-^/k^ + k^^ described by a 
Unearly polarized plane wave, 



Aph— e''*'-' sin (k • r+k^z — ojt) 



(5) 



whereas the DC applied field is Ae2;t= (0, — ci?t). For 
regularization we make use of a finite box L x L x Lz, 
so that momenta are discrete and the single photon's 
electric field is Eq = hoj/ (L^Lz). The unit vectors 



e'^^ = {— sin ip, cos ip) ; 



,(1) 



0; 



(6) 



e*-^^ = — cos6' (cos siniys) ; e^^-* 



smf 



describe polarizations that are conveniently chosen simi- 
larly to a recent calculation of electromagnetic emission 
due to thermal fluctuations 2l| . The vectors e*^^) and e*^^) 
represent the "in plane" and the "out of plane" polariza- 
tions, respectively. The electron and the hole wave func- 



tions are "^/fj^^*'' (p) V'n i^) ^-^d 
with spinors defined by 



V2L 



1 

-ie 



Z0 



; v{p') 



J4>' 



(7) 



with p+-^Aea;t=p (cos0, sin0) . 'ipni'^) wave func- 
tions of the confinement. The interaction with a pho- 
ton at time t happens when the momentum is minimally 
shifted due to the DC field. The Golden rule photon 
emission rate (for an "initial" electron with momentum 
p and a "final" hole p' and a photon of polarization A 
and momentum (k, k^)) is 



iVp(t)7V_p, it) (8) 



x5 (hvg [p + p') — huj) . 



In terms of Feynman diagrams of Quantum 
Electrodynamics [i^ it corresponds to the diagram 
in Fig. la. Here A'p (t) is the density of electrons in a 
certain momentum range produced by the electric field 



E, and A''_p/ (t) the density of holes (equal to that of 
the electrons at the opposite momentum due to particle 
- hole symmetry). The density calculated using the 
simple Landau - Zener creation rate expression for one 
of the flavours is [11 [11 [11 : 



iVp it) = e ipy) e [^Et - py) exp (-^P^) , (9) 

where G are the Heaviside functions. The transition am- 
plitude is given by 



(A) ^ 'iiv,ip+p')-u)j:W 



(10) 



dz e*'=^^^:(z)V'„,(z)<5(p + p'-k), 



where matrix elements J^^pi = i^P') ' e^^-^u (p) are 



T-(l) 



1 — cos ( 2(/3 — — 



(11) 



(2) 



2cos^ 



1 + cos ( — — 



Spectral emittance 

For tight confinement to the z = plane one should 
consider only the ground state n = n' = 0. Note that the 
perpendicular component of the wave vector k^ is " free" 
from conservation that prohibits the process in fully rel- 
ativistic QEd[23- The phase space for annihilation is 
very limited due to Vg << c, see Appendix A and leads 
to important simplifications. 
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Fig. 2. The spectral emittance in direction perpendicular 
to the graphene plane, k = 0. Polarizations are summed 

over. The omittance (in units of e^/t„;) for various 
frequencies (in units of t^^^ ) as function of ballistic time 
from O.ltni to lAtni- 
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Fig. 3. The emittance at various times (in units of 
as function of frequency. 

Let us define the spectral emittance per volume of the 
k - space (and area of the graphene flake) as 



M^(^) (p,p',k,fc„0 



dp 



(A) 



p,k-p 



(2^) 

X S{vg{p+\k-p\) -oj) 



^pA^k-p 



(12) 



where the integration over p' was performed using 
the delta function expressing the conservation of mo- 
mentum. We first study the frequency dependence 
of the radiation in the direction perpendicular to the 
graphene flake, k = 0. Multiplying with 4 for the spin 
and valley degeneracy, summing over the polarizations 

r p^p 



usmg uj 



~ 4 J , and integrating over p one obtains, 
cfcz, the spectral emittance 



(k = 0,a;,i) 



t/tn 



dp eitniuj/2+p) 



,exp[-2,r(t>V4-p^)] 



1/2 



(13) 



The spectral emittance, presented in Fig. 2 for various 
frequencies as function of time, increases linearly for t << 
a;"\ M{k = 0,uj,t) = ^^^g-'^tL'-'Va^ then rises sharply 

approaching a maximum at t = Lot^j2 and stabilizes at 



(k = 0,a;,t >> t,u) 



(14) 

where /q is the modified Bessel function. The asymp- 
totic value rises linearly with frequency, 7r~^a;, in the in- 
frared, reaches its maximum at w = tr^} and falls slightly 



to 



V2e 



in the ultraviolet. In Fig. 3 the emittance at 
various ballistic times is given as function of frequency. 



For each ballistic time the curve has two parts. The 
first follows the universal dependence given by Ea. (|T4l) . 
Therefore the frequency for observation of the Schwinger 
effect, not surprisingly, should exceed Wmin = y/^Evg/h^ 
that amounts to 3.6THz for E = lO^V/m, and bOTHz 
ioi E = 2 ■ lO^V/m. At a higher frequency cjniax ~ 2t/ 
t^i the emittance sharply drops. Therefore the frequency 
does not exceed 2tbai/t^i- 




Fig. 4a. The angle dependence of the intensity as 
function of the photon spherical angles 9 (half of the 
whole solid angle). Time is fixed at t = tni- The in - 
plane polarization. 




Fig. 4b. The out of plane polarization. 



Angular and polarization distribution 

Next we consider the angular and polarization depen- 
dence of the radiated power per unit area defined as the 
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spectral intensity, Ea (jl2l) . integrated over frequencies: 

£W (0, t) ^ r du'^M^^^ (k, -,t). (15) 

Performing integrations and simplifying, utilizing the 
small parameter v = Vg/c ~ 1/300 << 1, see Appendix 
B for details, one obtains: 

(16) 



are of the same order, while for 9 ^ 90° (close to the in 
- plane direction) the "out of plane" polarization, A = 2, 
dominates. On the other hand, for 0° << (p << 180° 
the picture is opposite. As expected, the unpolarized 
intensity. Fig. 2c, is less anisotropic; yet the radiation 
is somewhat depressed in the direction perpendicular to 
the current and close to the plane. 

The two-photon processes, see Fig. 2b, are suppressed 
by the factor of uqed = ~ Tjy compared to the one- 
photon process due to an additional vertex, while the 
phase space of two diagrams is of the same order; see 
Appendix C for details. 



£(2) (0,^,i) 



25/27r4c;i2 

■ sin^ (fi + —3- cos"' If 1 .(17) 



5tt 



The radiant flux from a flake of a /zm x fj,m size is 
4.7-10~'^^W, for E = W^V/m corresponding to the emis- 
sion rate of just 10 photons per second, yet for the high 
current samples^ with E = 2- Iff'V/m of the same area 
one gets a more significant output: the radiant flux is 
1.3- 10~^^Ty, corresponding to the emission rate of 3 • lO"' 
photons per second. 




Fig.4c Unpolarized light. 

The two quantities C'^^''^^ and their sum are presented 
for t = tni in the spherical plots Fig.4a-c, respectively. 
The radiated power is maximal in direction perpendicular 
to the graphene plane. For directions close to the azimuth 
angles (p ^ 0° and 180° (perpendicular to the electric field 
or current) at small polar angles 9 ( perpendicular to the 
graphene plane) the intensities of the two polarizations 



DISCUSSION 

Now we elaborate on a number of related issues and 
comment on the experimental feasibility of exploring the 
Schwinger phase physics. We start with a qualitative 
discussion of the rather unusual physics arising at strong 
applied fields, when the pair recombination becomes an 
important relaxation channel. 



Coulomb interaction and the formation of the 
neutral electron-hole plasma 

Even if the ballistic time and the relaxation time are 
very large, Schwinger's dynamical pair creation phase 
cannot persist for a long time at large field since density 
of charges (or both signs) becomes large. When the den- 
sity of quasiparticles reaches the order of Pp = lO^^cm"^ 
a neutral electron - hole plasma is created (like in 
some semiconductor systems under irradiation). In this 
state electrostatic interactions (despite being screened at 
large distances) become dominant, as was discussed ex- 
tensively in connection with electron - positron pairs cre- 
ation in Quantum Electrodynamics [20|. When electrons 
and holes are close enough they strongly attract each 
other effectively facilitating the recombination process. 
The rate therefore far exceeds the one calculated within 
perturbation theory in Section III. Let us first estimate 
when this state is achieved at experimentally accessible 
situations. 

Assuming the Schwinger pair creation rate, Eq.Q, the 
density will approach Pp at times of order tp oc p JE?!"^ . 



With a moderate field value oi E = 10*V/m the 



'plasma time' tp = 140ps ~ 400 t„/ exceeds the ballistic 
time of the L = lfj.m long sample (and probably also 
the relaxation time in current experiments on graphene). 
Yet with higher achievable fields E ^2- lO^V/m [l], the 
plasma time is reduced to tp = 40/s ~ 2 t„/ <ti thai- 
Therefore one expects that the " radiation friction" dissi- 
pation channel opens up: electron - hole pairs annihilate 
emitting photons, which take energy out of the graphene 
sheet and thus a new Ohmic behavior is reached. This 
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is roughly the bahistic time range for which the emis- 
sion was calculated in Section III. Of course, due to 
the Coulomb attraction enhancement of the recombina- 
tion the intensity becomes grossly underestimated in the 
plasma regime. The pair density will have to be re- 
calculated via Boltzmann equations; this will be done in 
a separate publication. One however might expect qual- 
itatively that this conductivity in the "second" Ohmic 
regime certainly exceeds and is likely to reach several 
times (72 ■ Data presented in ref.Q for clean samples, 
see Figs. 3 and 12 therein, is consistent with the linear 
(Ohmic) I-V curve at such conductivity value. However 
the experimental situation in the transport experiment 
is rather complex as discussed next. 



Experimental evidence of the pair creation in 
samples of mesoscopic size 

In this subsection we use the above " radiation friction" 
scenario to discuss whether there is a clear and unam- 
biguous signature of the Schwinger's pair creation phase 
in transport experiments. In a series of remarkable exper- 
iments the nonlinear I-V were measured at high electric 
fields of the order 2 • IQ^V/m at room temperature^. 
The samples were treated in such a way that, despite 
being non-suspended, the typical charge asymmetry did 
not appear and the Dirac point was accessed convinc- 
ingly at zero gate voltage (this demonstrates high qual- 
ity and is in variance with most samples on substrate). 
The I-V curves were studied in various high and low mo- 
bility samples (up to ^ = 7000cto^/ [Vs)) and effects of 
disorder were partially controlled by irradiating the sam- 
ples. The samples were L ~ 1 — 2fim long and rather 
narrow {W ~ 0.5/im) and the four-probe technique was 
applied. Although a nonlinear I-V dependence / cc V" 
with exponent a = 1.3— 1.5 was observed at Dirac point, 
surprisingly the nonlinearity disappeared in the highest 
mobility samples. In these experiments the ballistic time 
thai = '2ps is much larger than tni = 20/s at the high- 
est applied voltage of AV. Unfortunately the crossover 
voltages, namely when tni (V) = thai or 

Vni = ^ = 0.32my (18) 

were not probed since at room temperature fc^T = 
25meV. As argued in the previous subsection the ra- 
diation friction causes a second Ohmic regime and the 
I-V curve in the clean samples is expected to be linear. 
It is disorder that might have caused the observed non- 
linearity in irradiated samples. This requires an addi- 
tional theoretical study that includes the effect of pair 
recombination. As argued above it becomes as important 
as the Landau-Zcncr-Schwinger pair creation process at 
such currents. 



CONCLUSIONS 

To summarize, electron - hole pairs are copiously cre- 
ated via Landau-Zener-Schwinger mechanism near the 
Dirac points in graphene or similar 2D electronic systems 
by an applied electric field, provided the available ballis- 
tic time exceeds i„/, Eq.([T]). The recombination into pho- 
tons produces a characteristic signal proportional to E^^'^ 
at frequencies of order t~j that enables unintrusive and 
unambiguous experimental observation of the Schwinger 
phenomenon. The angular and polarization dependence 
of the emitted photons with respect to the graphene sheet 
was calculated. At very high currents and sufficiently 
long ballistic times the recombination process becomes 
greatly enhanced by the electron - hole attraction and 
the radiation becomes an elfective channel of dissipation, 
the radiation friction. 

The calculation can be trivially extended to any sys- 
tem with a Dirac point - like spectrum as double layer 
graphene and the recently synthesized family of materials 
called " topological insulators" in which surface exci- 
tations are similar to those in graphene with the notable 
exception of chirality. Schwinger's mechanism is also ex- 
pected in these materials since the mechanism does not 
involve chirality (left and right movers contribute equally 
to the emission rate of graphene). These materials have 
an advantage of not being strictly two dimensional, al- 
though ballistic times might be shorter at present. The 
driving current should not necessarily be DC, a sufficient 
condition is uJext « tni ■ Detectors of light (photon 
counters) in the microwave-infrared which are sensitive 
enough have recently been developed [1^. Hopefully 
Schwinger's pair creation rate formula can be directly 
tested using novel condensed matter materials endowed 
with relativistic fermion spectra. 
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APPENDIX A. PHASE SPACE OF THE 
ONE-PHTON PROCESS 

Planar electrons and holes are described by their mo- 
menta p and p' in the x — y plane, (see Fig. la and sub- 
section IIA for notations), while the momentum of the 
photon k is three dimensional. The conservation of the 
in-plane momentum and the conservation of energy read: 
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p + p'=k; Vg{p + p') = ck. 



(19) 



The momentum in the z direction is not conserved; it is 
balanced by the elasticity of the graphene flake. Since 

k ~ \J k1 + |k|^ > |k| , in terms of p and p' one has the 
inequality 



^''(p+p')'-|p + p'r >o, (20) 

or, in the polar coordinates, 

(1 - v^) (p^ +p'^) + 2 [cos (0 - 0') - w^] pp' < 0. (21) 

With the above constraint the condition that p' has real 
solutions leads to 



cos I 



-0') < -l + 2i;2. 



(22) 



As t; ^ 1, it is obvious that 0' — is very close to tt. By 
defining = tt — (^' — 0) , one see the above condition 
simplifies to 



2v < A(f>< 2v. 



(23) 



Substituting the above result back into Ea. ((2T|) it can 
be seen that p' is very close to p. By introducing Ar — 
1 — p' /p, the condition becomes 

Ar^ + A^^ <4z;2. (24) 

Therefore the allowed region is a disk of radius 2v. 



(1) 



4 [cos {(j) — if) — V sin t 



(2) 
P,P' 



(l — 2w sin Q cos (i/i — (/j) + sin^ B) 
4cos^ 0sin^ (0 - ip) 



(1 — 2usin0cos ( 



V?) + 1-2 sin^ B) 



.(26) 



Moreover, there is an additional factor from the delta 
function: 



J(0) 



{\ — 2v sin Q cos (0 — (/s) + sin^ ( 
2vg (1 — 2z;sin0cos (</> — (^))^ 



(27) 



The Jacobian of the transition to polar coordinates is 



P(0) 



Lij {\ — sin^ I 



2vg (1 — 2v sin cos (0 — kp)) 



(28) 



In view of the step functions for the LZS density, Eq.(l9]), 
the following possibilities occur: 

a. into the forward direction (positive projection on the 
electric field) < <p < tt : 

The conditions imposed by the step functions are 
< —py and < ^t + Py ^ ky. In terms of the po- 
lar coordinates in the momentum space 

(i) for < < 1 + V sin 9 sin ip, the allowed regions 
are A_|_ — i^q < < or — tt < <j) < — tt + 4>q — A_|-, with 



arcsm 



2t 



A+ = i;sin6'tan(/)g (cos ((/3 + (/^o) + 2sin(^csc(/)oX29) 



APPENDIX B. DERIVATION OF THE 
AMPLITUDE AND SPECTRAL EMITTANCE 

The Golden rule photon emission rate (for an " initial" 
electron with momentum p, a "final" hole p' and a pho- 
ton of polarization A and momentum (k, k^)) is given by 
Eq.®. Correspondingly the rate defined in Eg. p^ is 



2 2 

{2^f 

9 9 

{2.t 



dp 



p,k p 



NpNi^^pS {vg (p + \k - p\) ~ Lo) 



p {(j)) d(i) 



P.k-p 



iVpiVk_p • J (</)) . 



(25) 



Because of the delta function, the integration over p 
imposes the condition of energy conservation. As a result, 
the squares of the matrix elements, Eq. (|lip simplify: 



(ii) for 1 -I- w sin 9 sin < -^^r- , the allowed region for (j) 
is — TT < (p <0. 

b. into the backward direction (negative projection on 
the electric field) —Tr<tp<0: 

The conditions imposed by the step functions now be- 
come < —py + ky and < ^t + Py. In terms of the 
polar coordinates in the momentum space 

(i) for -~2vsm9sm(p < < 1 — w sin 6* sin the 

allowed regions are — ((/)q — A_) < < 2t;sin0sin(^ or 



-TT — 2'y sin 9 sin Lp < (p < 



A_, where 



A_ = V sin 9 tan (p^ cos {ip + (j)^) . 



(ii) for 1 — w sin 9 sin Lp < 



(30) 



the allowed region for cj) 
is — TT — 2v sin 9 sin (p < (f) < —2v sin 9 sin ip. 

Since v = Vg/c ^ 1/300 << 1, one can neglect higher 
order correction in v. The conditions for case (i) simplify 



into < 



< 1 and 



< 



< or -TT < 



< 



8 



Therefore 



and 



M^'^ (6 



UJ 



87r4 



vtj ■ e 



,(A') 

{v(T ■ e'^) 



) cos^ ip 



Sir* 



cos 



) siii^ (f - 



sm 



)COS' 



' (fi) exp 



(fi^ exp 



■ {k-p')-a--(k~p') 

■ (p-fc)-iT-(p-k) 

v[{p-kf-{p-kf\ 



vcr ■ e 



,(A') 



2.2 



2 



cos 



Similarly, the conditions for case (ii) simplify: 
and -TT < 6 <0. Thus 



2t 



> 1 



The outcome of the integral in Eq. ([T^ is dictated by 
the size of the phase space and by the powers of the small 
parameters v = Vg/c and agsD = e^/ (he). The power 
, qfj i; is the same as in the one-photon case, see Eq. ([T0| . 
Yfit there appears an additional power of aqED- 

Let us estimate the phase space. The conservation of 
momentum and energy for the two-photon process imply: 



k + k'; 



.{p + p')^c{k + k'). 



(35) 



167r4 

o2 



UJ dcf) cos^ {(j) — (p) exp 



-o<i^ COS^ 



k'y) one still has 
> 0. Just like in the 



V4 



327r3 

(e, 



lo 



cos 2ipli 



{k + k'f > {k^^k'^f + (A:y + 

the inequality v'^{p + p'Y — |p + p'|^ 
one-photon case (Appendix A), utilizing v « I leads 
to the constraints — 0' « tt and p' « p. The phase 
i space of the two-photon case is thus roughly of the same 
order as that of the one-photon process. Therefore this 
two-photon process can be neglected. 



167r4 
327r3 



we 



V4, 



(ys) exp 



o;^ cos^ I 



cos 2ipli 



These expressions lead to the final results for the spectral 
emittance Eq.(ini) and Eq. p^ . 

Next we consider the angular and polarization de- 
pendence of the radiated power per unit area defined 
as the spectral intensity {9,(f,t), Eq.([T5]). In the 
above formulae one can first integrate over uj in the range 
< UJ < -3^ CSC (p, and then integrate over (p in the region 

'■nl 

[O, f ] . This leads to Eqs.dH]) and ^ for the luminosity 
integrated over frequencies. 



APPENDIX C. PHASE SPACE OF THE 
TWO-PHOTON PROCESS 

As in QED [^O], the two-photon diagram. Fig. lb, 
gives rise to the following S - matrix element 



,(A,V) 



le 



{p,p',k,k'] 

(A,A') 



A^ujuj'ee' 



(33) 

{2t:)U {vg{p + p')-c{k + k')), 



where 



(A,A') 



(-p')s(^'^')„(p) 



(34) 
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